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Abstract. A system of singularly perturbed ordinary differential equations of first order with 
given initial conditions is considered. The leading term of each equation is multiplied by a small 
positive parameter. These parameters are assumed to be distinct and they determine the different 
scales in the solution to this problem. A Shishkin piecewise-uniform mesh is constructed, which is 
used, in conjunction with a classical finite difference discretization, to form a new numerical method 
for solving this problem. It is proved that the numerical approximations obtained from this method 
are essentially first order convergent uniformly in all of the parameters. 

Key words, linear dynamical system, multiscale, initial value problem, singularly perturbed, 
finite difference method, parameter-uniform convergence 

AMS subject classifications. 65L05, 65LI2, 65L20, 65L70 

1. Introduction. We consider the initial value problem for the singularly per- 
turbed system of linear first order differential equations 

(1.1) Eu'{t) + A(t)u(t) = f(t), te(0,T], u(0) given. 

Here u is a column n- vector, E and A(t) are n x n matrices, E = diag(e), e = 
(e\, . . . , £„) with < Ei < 1 for all i = 1 . . .n. For convenience we assume the 
ordering 

£i< ... <e n . 

These n distinct parameters determine the n distinct scales in this multiscale problem. 
Cases with some of the parameters coincident are not considered here. We write the 
problem in the operator form 

Lu = f, u(0) given, 
where the operator L is defined by 

L = ED + A(t) and D = — . 
We assume that, for all t E [0, T], the components a,ij(i) of A(t) satisfy the inequalities 

n 

(1.2) a u (t) > y2\aij(t)\ for 1 < i < n, and a^t) < for i ^ j. 

.1 = 1 

We take a to be any number such that 

n 

(1.3) < a < min (Va,,(()). 

te(o,i] 

l<i<« 3=1 
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We also assume that T > 2max ! (£ i )/a, which ensures that the solution domain 
contains all of the layers. This condition is fulfilled if, for example, T > 2/ a. We 
introduce the norms || V ||= maxi<fc< n T4| for any n-vector V, \\ y \\— sup 0<(<T \y(t)\ 
for any scalar- valued function y and || y \\ = maxi<fc<„ || y^ || for any vector- valued 
function y. Throughout the paper C denotes a generic positive constant, which is 
independent of t and of all singular perturbation and discretization parameters. Fur- 
thermore, inequalities between vectors are understood in the componentwise sense. 

The plan of the paper is as follows. In the next section both standard and novel 
bounds on the smooth and singular components of the exact solution are obtained. 
The sharp estimates in Lemma 2.4 are proved by mathematical induction, while an 
interesting ordering of the points tij is established in Lemma 2.6. In Section 3 the 
appropriate piecewise-uniform Shishkin meshes are introduced, the discrete problem 
is defined and the discrete maximum principle and discrete stability properties are 
established. In Section 4 an expression for the local truncation error is found and 
two distinct standard estimates are stated. In the final section parameter-uniform 
estimates for the local truncation error of the smooth and singular components are 
obtained in a sequence of lemmas. The section culminates with the statement and 
proof of the parameter-uniform error estimate, which is the main result of the paper. 

The initial value problems considered here arise in many areas of applied mathe- 
matics; see for example pQ. Parameter uniform numerical methods for simpler prob- 
lems of this kind, when all the singular perturbation parameters are equal, were 
considered in [4]. A special case of the present problem with n — 3 was considered 
in [3J, which also contains numerical results confirming the theory. For this reason 
further numerical validation is considered to be unnecessary. A general introduction 
to parameter uniform numerical methods is given in [2] and [7]. 

2. Analytical results. The operator L satisfies the following maximum princi- 
ple 

Lemma 2.1. Let A(t) satisfy il.ty) and il.S\) . Let ip(t) be any function in the 
domain of L such that tp(0) > 0. Then Lip(t) > for all t G (0, T] implies that 
> for all t G [0, T] . 

Proof. Let i*, t* be such that ipi*(t*) = min^t tpi(t) and assume that the lemma is 
false. Then ^i* (t*) < . From the hypotheses we have t* ^ and ip[* (i*) < 0. Thus 

n 

(Lip(t*))i* = Ei-Vv (**) + a i',i* (t*)ipi* (t*) + a **An^(n 

j=l, j^i* 

n 

which contradicts the assumption and proves the result for L. □ 

Let A(t) be any principal sub-matrix of A(t) and L the corresponding operator. 

To see that any L satisfies the same maximum principle as L, it suffices to observe 
that the elements of A(t) satisfy a fortiori the same inequalities as those of A(t). 

We remark that the maximum principle is not necessary for the results that follow, 
but it is a convenient tool in their proof. 
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Lemma 2.2. Let A(t) satisfy ll.Sfy and il.3\) . If ip(t) is any function in the 
domain of L then 

|| m || < max 1 1| $0) ||, i || U ||| , t G [0,T] 



Proof. Define the two functions 

± W=max{||^(O)||,i||^||}e±^), 
a 

where e = (1, ... ,1)' is the unit column vector. Using the properties of A it is 
not hard to verify that 0^(0) > and L0 ± (i) > 0. It follows from Lemma T2.ll that 
&±{t) > for all t G [0,T]. □ 

The Shishkin decomposition of the solution -u of (1) is given by u = v + w where 
v is the solution of Lv = f on (0, T] with v(0) = ^4 _1 (0)/(0) and w is the solution 
of Lw = on (0, T] with w(0) = u(0) — v(0). Here v , w are, respectively, the smooth 
and singular components of u . 

The smooth component v of u and its derivatives are estimated the following lemma, 
which gives bounds showing the explicit dependence on the inhomogeneous term and 
the initial condition. 

Lemma 2.3. Let A(t) satisfy and il.3\) . Then there exists a constant C, 

independent o/£,-u(0) and f , such that 

||*||<<7||/||, IKII<c(ll/ll + llfll) 

and, for all 1 < i < n, 

II e iV 'l ||<C(||/|| + ||/'||) 



Proof. We introduce the two functions i/' ± (^) = C||/l|fT± v(t) where e is the unit 
column vector. Noting that v(0) = A _1 (0)/(0), it is not hard to see that V' ± (0) > 
and Lif ± (t) > 0. It follows from Lemma |2~T1 that ^(t) > for all t e [0, T] and 
so || v || < C || / || . To estimate the derivative we now define the two functions 
<£±(t) = C(\\f\\ + \\f'\\)e±v'{t). Since v'(0) = and Lv' = /' - A'v, it may be 
verified that ^(0) > and L$ ± (t) > 0. Again by Lemma ETC! we have $ ± (t) > 0, 
which proves the result. Finally, differentiating the equation £iv[ + {Av)i = fi and 
using the estimates of v and v' , we obtain the required bound on Eiv" □ 

We define the layer functions Bi, 1 < i < n, associated with the solution u by 

Bi(t) = e- at/e % t £ [0,oo). 

The following elementary properties of these layer functions, for all 1 < i < j < n, 
should be noted: 

(i) Bi(t) < Bj(t), for all t > 0. 

(ii) Bi(s) > Bi(t), for all < s < t < oo. 

(iii) B % (0) = 1 and < B t {t) < 1 for all t > 0. 

Bounds on the singular component w of u and its derivatives are contained in 



4 



S. Valarmathi and John J.H. Miller 



Lemma 2.4. Let A(t) satisfy ll.Sfy and il.3\) . Then there exists a constant C, such 
that, for each t € [0, T] and i = 1, . . . , n, 

K(t)| < cB n (t), K(*)| < c^^S, \ £i w'i(t)\ < cj^^r 1 - 

q—l ^ q— 1 ^ 

Proof. First we obtain the bound on w. We define the two functions ip^ = 
CB n e± w. Then clearly ^±(0) > and bip ± = CL(B n e). Then, for i = 1, . . . , n, 
(Lijj ± ) l = C{J2] =1 a^j - af-)B„ > 0. By Lemma HO # > 0, which leads to the 
required bound on w. 

To establish the bound on w' we begin with the n th equation in Lw = 0, namely 

e n w' n + a nA w 1 + . . . + a n>n w n = 0, 

from which the bound for i = n follows. We now bound w[ for 1 < i < n — 1. We 
define p = (wi, . . . , w n -i) and, taking the first n — 1 equations satisfied by w, we get 

Ap = g, 

where A is the matrix obtained from A by deleting the last row and column and the 
components of g are gk — —CLk,nW n for 1 < k < n — 1. Using the bounds already 
obtained for u> we see that g is bounded by CB n (t) and its derivative by C B ™^ . 
The initial condition for p is p(0) = u(Q) — u°(0), where if is the solution of the 
reduced problem if = A^ 1 f ', and is therefore bounded by C(|| u(0) || + || /(0) ||). 
Decomposing p into smooth and singular components we get 

p = q + r, ff = <f + ? . 

Applying Lemma 12.31 to q, from the bounds on the inhomogeneous term g and its 
derivative <f, we conclude that || (fit) ||< C B ^ ^ . We now use mathematical induc- 
tion. We assume that Lemma [2^41 is valid for all systems with n — 1 equations. Then 
Lemma |2 .41 applies to r and so, for i = 1, . . . ,n — 1, 



En-l 



Combining the bounds for qi and we obtain 



\ P m\<c(^+...+^). 



Recalling the definition of p this is the same as 



We have thus proved that Lemma 12.41 holds for our system with n equations. Since 
Lemma 12.41 is true for a system with one equation, we conclude by mathematical 
induction that it is true for any system of n > 1 equations. 

Finally, to estimate the second derivative, we differentiate the i equation of the 
system Lw — to get 

e k w'l = -(Aw- + A'w)i 
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and we see that the bound on w'! follows easily from the bounds on w and w' . □ 
Definition 2.5. For each 1 < i =/= j < n we define the point tij by 

( 2 j\ B i(tj,j) _ Bjjtjj) 

Si Sj 

In the next lemma it is shown that these points exist, are uniquely defined and have 
an interesting ordering. Sufficient conditions for them to lie in the domain [0, T] are 
also provided. 

Lemma 2.6. For all i,j with 1 < i < j • < n the points tij exist, are uniquely 
defined and satisfy the following inequalities 

(2.2) e- 1 B l {t)>e] 1 B ] (t) te[0,t«) 
and 

(2.3) el l Bi{t) < eJ l Bj{t) t e (ty.oo). 
In addition the following ordering holds 

(2.4) tij < tj+ij, if i + 1 < j and < if i < j 
and 

(2.5) <£j/2 implies that £ (0,T] for all i < j. 



Proof. Existence, uniqueness, (|2.2[) and (|2.3[) all follow from the observation that 
for i < j we have e, < £j and the ratio of the two sides of (|2.ip , namely 

Bi(t) £o Sj , .1 1 

1 -o k = fexp (-<---)), 

is monotonically decreasing from > 1 to as i increases from to oo. 
Rearranging (|2.ip gives 

Mi)-Mi) 

v Si ej I 

Writing Ek = exp(— pk) for some pk > and all k gives 

t _ - Pj 

1,1 a(exppi - exppj) 

The inequality t% j < it+i.j is equivalent to 

Pi - Pj Pi+i - Pj 



euppi — exppj expp, + i — exppj 
which can be written in the form 

(Pi+i - Pj) exp(p l -pj) + {pi - p i+ x) - (pi - p^ exp(p J+ i ~Pj)> 0. 
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With a = pi — pj and b = Pi+\ — Pj it is not hard to see that a > b > and 
a — b = pi — Pi+\. Moreover, the previous inequality is then equivalent to 

exp a — 1 exp 6—1 
a b 

which is true because a > b and proves the first part of (|2.4| . The second part is 
proved by a similar argument. 



Finally, to prove (|2.5[) it suffices to rearrange (|2.ip in the form 



t 



h3 a(±-±Y 

Since T > § and Si < % it follows that ln(J-) < |f and e (0, T]. □ 

3. The discrete problem. We construct a piecewise uniform mesh with TV 
mesh-intervals and mesh-points {ti}^L by dividing the interval [0, T] into n + 1 sub- 
intervals as follows 

[0,T] = [0,a 1 ]U(a 1 ,a2]U...(a n -i,(T n ]l>((T n ,T] 

Then, on the sub-interval [0, 01], a uniform mesh with ^ mesh-intervals is placed, 
and similarly on (<7j, crj+i], l<i<n — 1, a uniform mesh with 2 „^ +1 mesh-intervals 
and on (cr n , T] a uniform mesh with -y mesh- intervals. In practice it is convenient to 
take N = 2 n k where k is some positive power of 2. The n transition points between 
the uniform meshes are defined by 

<Ji = mm{— — , — In TV} 
2 a 

for i = 1 , . . . , n — 1 and 

T e n 

<j„ = min{ — , — In TV}. 



Clearly 



T 

< cri < ... < a n < -. 



This construction leads to a class of 2 n piecewise uniform Shishkin meshes Mr, where 



6' 

6 denotes an n-vector with bi = if Ui = and 6^ = 1 otherwise. Writing 

<5j = tj — we remark that, on any Mr, we have 

(3.1) Sj < CTV -1 , l<j<TV 
and 

(3.2) < Ce; In TV, 1 < i < n. 

On any Mr we now consider the discrete solutions defined by the backward Euler 
finite difference scheme 



ED-U + A(t)U = f, U(0) = u(0), 
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or in operator form 

L N U — f, U(0)=u(0), 

where 

L N = ED~ + A(t) 
and D~ is the backward difference operator 

"3 

We have the following discrete maximum principle analogous to the continuous case. 

Lemma 3.1. Let A(t) satisfy kl.3!\) and \ 1.3\) . Then, for any mesh function ^ , the 
inequalities *(0) > and L N *(tj) > for 1 < j < N, imply that $(tj) > 
/orO < j < N. 

Proof. Let be such that Vi*(ij») = min^- Vi(tj) and assume that the lemma 
is false. Then Vi*(tj*) < . From the hypotheses we have j* ^ and Vi*(tj») — 
Vi*{tj*-x) < 0. Thus 

(L*?(fj.))i- = e .. ^fe-)-Ji-fe--i) + a-.^.^.)^.^..) + 

n 

<^i'(*j*) I] a • »• 

fc=l fc^i* 

which contradicts the assumption, as required. □ 

An immediate consequence of this is the following discrete stability result. 
Lemma 3.2. Let A(t) satisfy and il.3\) . Then, for any mesh function ^, 

II II < max 1 1| %) ||, 1 || L N * ||| , < j < N 

Proof. Define the two functions 

9 ± (t) = max{||*(0)||, i||i>*||}e± ¥(*) 
a 

where e = (1, ... , 1) is the unit vector. Using the properties of A it is not hard to 
verify that 6^0) > and L^9 ± (^) > 0. ft follows from LemmaOthat 6 ± (^) > 
for all < j < N. Q 

4. The local truncation error. From Lemma 13.21 we see that in order to 
bound the error || U — u || it suffices to bound L N (U — u). But this expression satisfies 

L N (U - tf) = L N (U) - L N (u) = f- L N (u) = L(u) - L N \u) 

= (L — L N )u = -E{D~ - D)u, 
which is the local truncation of the first derivative. We have 

E(D- -D)u = E(p- -D)v + E{D~ - D)w 
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and so, by the triangle inequality, 

(4.1) || L N (U - u) \\<\\ E(D~ - D)v || + || E(D~ - D)w \\ . 

Thus, we can treat the smooth and singular components of the local truncation error 
separately. In view of this we note that, for any smooth function ip, we have the 
following two distinct estimates of the local truncation error of its first derivative 

(4.2) \(D~ - D)tp(tj)\ < 2max|V''(s)| 

seij 

and 

(4.3) |(D--D)^)|<|ma X |^"( S )|, 

where Ij — [tj-i,tj]. 

5. Error estimate. We now establish the error estimate by generalizing the ap- 
proach based on Shishkin decompositons used in [5] . For a reaction-diffusion boundary 
value problem in the special case n = 2 a parameter uniform numerical method was 
analyzed in [6] by a similar technique and in the general case in [5] using discrete 
Green's functions. 

We estimate the smooth component of the local truncation error in the following 
lemma. 

Lemma 5.1. Let A{t) satisfy \1.2\) and fl.3\) . Then, for each i = 1, . . . , n and 
j = 1, . . . , N, we have 

\si{D- - D) Vi (tj)\ <CN-\ 

Proof. Using (|4.3[) . Lemma |2~31 and (|3.ip we obtain 

\£i(D- — D)vi(tj)\ < C^max|£X'( s )l < Cfy < CN^ 1 

seij 

as required. □ 

For the singular component we obtain a similar estimate, but in the proof we must 
distinguish between the different types of mesh. We need the following preliminary 
lemmas. 

Lemma 5.2. Let A{t) satisfy M.2\) and \1.S\) . Then, for each i = 1, . . . , n and 
j = 1, . . . , N, on each mesh Mg, we have the estimate 

\e i (D--D)w i (t j )\<C^. 

El 

Proof. From (|4.3j) and Lemma E2J we have 

\ Si (D- -D)wi{tj)\ < 
< 
< 



C5j maxs^. \e i w i (s)| 
l^ q =i — T q — 

£l 



as required. □ 
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In what follows we make use of second degree polynomials of the form 

E(t - t e ) k (fe) . . 
jfci w i '(*«). 

A:=0 

where 6 denotes a pair of integers separated by a comma. 

Lemma 5.3. Let A(t) satisfy fl.tfy and fl.3\) . Then, for each i = 1, . . . , n, 
j = 1, . . . , N and k = 1, . . . , n — 1, on eac/i mes/i Mg with b k = 1, there exists a 
decomposition 

fc+i 

771—1 

/or which we have the following estimates for each m, 1 < m < k, 
\sM, m (t)\ < CB m (t), \e lW 'l m (t)\ < 

Cm 

and 

B q (t) 



\e lW 'l k+1 (t)\ < C 

Furthermore 



q =k+i £ i 



\ei(D- - D)w i (t j )\ < C{B k {t^ x ) + -L-). 

£k+l 

Proof. Since b k = 1 we have e k < £fc+i/2, so t k .k+i € (0, T] and we can define the 
decomposition 

fe+i 

= ^ w i,m, 

where the components of the decomposition are defined by 



and for each m, k > m > 2, 



Pi;k,k+1 On [0,tfc,fc + l' 

u>,- otherwise 



ft;m-l.m On [0, t m „i !ln ) 

w * ~ L g=m +i otherwise 
and 

k+l 

m,i — m - ^2 w i,q on [°i T l- 

From the above definitions we note that for each m, 1 < m < k, Wi -m = on [t m)m +i, T]. 
To establish the bounds on the second derivatives we observe that: 
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in [tk.k+i, T], using Lemma l2"T4l and t > t^k+x, we obtain 

| £ ,< fe+1 (t)l = N<(i)i<cf:^<c £ 

g=l « «=fc+l 9 

in [0, tfc k+i], using Lemma 12.41 and t < tk t k+x , we obtain 

_„// (+\\ - \ e ,JI (4 \\<-X^ ^gfa.fc+l) ^ -Bq^MH-l) -Bgft) 

£ ^,fc+lWI - Fi^i l*fc,fc+Ul < 2^ £ - £ - _ £ 

q=l q q=k+X q q=k+X q 

and for each m = k, ... , 2, we see that 
in [t m . m +i,T}, w" m = 0; 

in [t m _i, m) im,m+i]) using Lemma T2. 41 we obtain 

ki< m (*)l < + £ | £i <,(t)| <cJ2^-< 

g— m+1 

in [0, i m _i, m ], using Lemma HOI and i < t m -x, m , we obtain 

I \ i ^- ^V 1 Bq{t m _x,m) , f~(Bm(t m -Xm) , r ,B m (t) 
(t)\ = \EiWi {t m -x,m)\ < C < 6 < C ; 

9=1 



in[t 1)2) T], < x =0; 

in [Ojti^], using Lemma HOI 



< W(*)i + E i e K«(*)i <^£^< c^. 

n ^ Eq £\ 

q=2 9=1 H 



For the bounds on the first derivatives we observe that for each to, 1 < to < k 
in [t m , m +x,T], w' lm = 0; 

in [0, t m<m+ x] J* m - m+1 ei w'l m (s)ds = eiw' i m (t mjn+1 ) - Siw' i m (t) = -eiw' i m (t) 
and so 

|£i< m (i)|< / ' + \e t w'l m (s)\ds < — f ' + B m (s)ds<CB m (t). 
Finally, since 

k 

\£i{D- - D)w i {t j )\ < J2 \ £ *( D ~ - £>Km(*j)| + h(D- - D)w itk+1 (tj)\, 

m=X 

using (|4.3[) on the last term and (|4.2p on all other terms on the right hand side, v 
obtain 

k 

\si(D- -D)wi(tj)\ < C(Y^ max|<?iW- )m (s)| + ^- max | £i u;^ +1 (s)|). 
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The desired result follows by applying the bounds on the derivatives in the first part 
of this lemma. □ 

Lemma 5.4. Let A(t) satisfy U.2\) and U.3\) . Then, for each i = 1, . . . , n and 
3 = 1, . . • , N, on each mesh Mg, we have the estimate 

\ £i (D- - D)w i (t j )\ < CBnitj-i). 

Proof. From (|4.2|) and Lemma \2. 41 for each i = 1, . . . , n and j = 1, . . . , N, we 
have 

|£i:(-D~ - D)wi(tj)\ < Cmax se /j£i^(s)| 
< CB n (tj-i) 

as required. □ 

Using the above preliminary lemmas on appropriate subintervals we obtain the 
desired estimate of the singular component of the local truncation error in the follow- 
ing lemma. 

Lemma 5.5. Let A(t) satisfy il.ty) and 11. 3\ ). Then, for each i = 1, . . . , n and 
3 = 1, • ■ ■ , N , we have the estimate 

\ei(D- - D)wi(tj)\ < CN- 1 IniV. 

Proof. We consider each subinterval separately. 
In the subinterval (0, a{\ we have Sj < CN^ai. On any mesh Mg, using Lemma HT2"! 
we get \si(D- - D)wi(tj)\ < CA^ 1 ^- < CN^laN. In the subinterval (a 1 ,a 2 \ we 
have 5j < CN~ 1 a 2 . 

On any mesh M-g with 6i = 0, we have o~2 = 2a%. Using Lemma 15.21 we get 
\ei(D~ -D)wi(tj)\ < CN' 1 ^ < CN- 1 ^ < CN^ItlN. 

On any mesh Mj with b\ = 1, we have <j\ — — IniV. Using Lemma [5.31 with 
k = 1 we get \ei{D~ - D)w i (t j )\ < C(B 1 (a 1 ) + A^ 1 ^) < CAT 1 IniV. 
In a general subinterval (<7 m , <7 m +i] 2 < m < n — 1, we have Sj < CN~ cr m +i- 

On any mesh with £> 9 = 0, g = 1, ... , m, we have cj m +i = Ccri. Using 
LemmaOwe get \si(D~ - D)wi{tj)\ < CN- 1 ^^- < CN' 1 ^ < CN^lnN. 

On any mesh Mg with b\ = 1, b q = 0, 9 = 2, ... , m, we have <Ti = — In iV, cr m +i = 
Ccr 2 . Using Lemma 15.31 with k — 1 we get 1^(1?^ — Z?)uii(tj)| < C{B 1 (a m ) + 
N- 1 ^^) < C(Bi((7i) +iV- 1 fJ) < CAT" 1 In AT. 

On any mesh Mr with 6^ = 1, b q = 0, g = fc + 1, ... , m, we have <jfc = 
^MnAT, a rn+ i = C<Jk+i- Using Lemma l5~3"l with general k we get \ej,{D~— D)wi(tj)\ < 
C(B k (a m ) + A^ 1 ^- 3 -) < C(B k (a k ) + A^ -1 ^^-) < CN^ 1 \\iN. 

c k-\- 1 c fc+ 1 

On any mesh with b m — 1, we have cr m = ^InN. Using Lemma I5T41 we get 
\ Ei (D- - D) Wi (tj)\ < CN^BUam) < CN^ 1 In N° 
In the subinterval (er„, T] we have Sj < CN^ 1 . 

On any mesh with & g = 0, q = 1, ... ,n, we have 1/ei < ClnA^. Using 
LemmaHlwe get |e<(D- - D) Wl {t )\ < CN^/si < CN^laN. 

On any mesh Mg with b\ = 1, 6 g = 0, q = 2, ... , n, we have <7i = — In A^, l/e2 < 
ClnA^. Using Lemma [Ol with k = 1 we get |e 4 (L>~ - D)^^)! < C(Bi(cr n ) + 
/e 2 ) < C(Bi(cti) + N- 1 /e 2 ) < CN' 1 InN. 
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On any mesh Mg with 6^ = 1, b q = 0, q = k + 1, ... , n, 2 < fc < n — 1, 
we have = — In TV, 1/ek+i < ClnJV. Using Lemma [5.31 with general k we get 
\e l {D--D)w l {t a j )\ < C( J B fc (a„)+^- 1 /e fc+1 ) < C(B fe (a fe )+iV- 1 /e fc+1 ) < CN^laN. 

On any mesh with 6„ = 1, we have cr n = — IniV. Using Lemma T5.4I we get 
|ei(D - - D)wi(tj)\ < CN- l B n {a n ) < CW^ 1 InN. " 

It is not hard to verify that on each of the n + 1 subintervals we have obtained 
the required estimate for all of the 2™ possible meshes. □ 

Let u denote the exact solution of and U the discrete solution. Then, the 
main result of this paper is the following e-uniform error estimate 

Theorem 5.6. Let A(t) satisfy tl.ty) and fl.3\) . Then there exists a constant C 
such that 

|| U-u \\< CN^hxN, 

for allN > 1 

Proof. This follows immediately by applying Lemmas 15.11 and 15.51 to (|4.1[) and 
using Lemma l3~2l □ 
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